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Abstract 



Two standard models of sol-gel transition are revisited here from the point of view 
of their fluctutations in various moments of both the mass-distribution and the gel- 
mass. Bond-percolation model is an at-equilibrium system and undergoes a static 
second-order phase transition, while Monte-Carlo Smoluchowski model is an off- 
equilibrium one and shows a dynamical critical phenomenon. We show that the 
macroscopic quantities can be splitted into the three classes with different scaling 
properties of their fluctuations, depending on wheather they correspond to : (i) non- 
critical quantities, (ii) critical quantities or to (iii) an order parameter. All these 
three scaling properties correspond to a sing le form : < M > 5 P(M) = $((M- < 
M >)/ < M > s ) , with the values of 5 respectively : = 1/2 (regime (i)), 7^ 1/2 and 
1 (regime (ii)), and = 1 (regime (iii)). These new scalings are very robust and, in 
particular, they do not depend on the precise form of an Hamiltonian. 

PACS numbers: 05.40.+j,82.20.-w,82.20.Mj 
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I. REVERSIBLE AND IRREVERSIBLE AGGREGATION MODELS 



Below we deal with a system of monomers (the basic units) which can aggregate to form 
the connected clusters. A monomer is considered as a cluster of mass 1 (this is the mass 
unit). Moreover, monodispersity of the monomeric mass as well as the mass conservation 
during the aggregation is assumed, i.e., the cluster-masses are always integer numbers in 
this approach. 

The basic sol-gel transition is the appearance at a finite time of an infinite cluster, called 
the gel. 'Infinite' means here that a finite fraction of the total mass of the system belongs to 
the gel. Note that this definition is applicable both to finite as well as to infinite systems, 
but contextually and historically, the tools used to study this behaviour have been defined 
somewhat differently in these two cases [T] . For finite systems, the following moments of the 
number-mass-distribution n s ( i.e. the number of clusters of mass s) are introduced : 

s 

where the summation is performed over all clusters except the largest one. The superscript 
' in Eq. ([]]) recalls this constraint on the allowed values of s . Consequently, the mass of 
gel-fraction is just iV — M[ with : 

N = J2 sn s , 

all s 

the total mass of the system. For infinite systems, the following normalized moments of 
the concentration-mass-distribution c s (i.e. the concentration of clusters of mass s) are 
introduced : 

m' k = J2s k c s , (2) 

s 

where the summation runs over all s . Generally, the concentrations are normalized in a 
special way : 

n 

c s = limAr^oo-^ , 

and not in a direct relation with the volume. Consequently, the probability for a monomer 
to belong to the gel-cluster is just equal to 1 — m' 1 . 

These different definitions of moments are equivalent in the sense that they are connected 
through the relation : 

, r K 

m k = hmjv^oo— . 

In this paper, we investigate scaling behaviours of distributions of mass-distribution- 
moments for the two 'classical' approaches to the sol-gel transition : the reversible model 
(the percolation M) and the irreversible model (the Smoluchowski equations In the 
former case, diffusion is unimportant and the reaction between clusters is the relevant step, 
whereas in the latter case, on the contrary, the reaction between clusters is unimportant 



2 



and the diffusion is relevant. In this sense, these two models are believed to belong to the 
two different classes, and all other models of aggregation are suspected to belong to one of 
them. For example, statistical Flory-Stockmayer theory |3| should behave as the percolation 
model for the scaling properties, since it corresponds to the calculation 'at equilibrium'. This 
'universality' in diverse sol-gel situations is at present a guess. 

Percolation model |IJ can be defined as follows : in a box (a part of the regular lattice), 
each site correponds to a monomer and a proportion p of active bonds is set randomly 
between sites. This results in a distribution of clusters defined as ensemble of sites connected 
by active bonds. For a definite value of p, say p cr , a giant cluster almost surely spans all the 
box. For example, in the thermodynamic limit when the size of box becomes infinite (this 
limit will be denoted by 'lim' in the following), a finite fraction of the total number of vertices 
belongs to this cluster. Therefore, we get the results : m! x = 1 for p < p cr and m! x < 1 for 
p > p cr . Moreover, m! x is a decreasing function of the occupation probability. This typical 
behaviour is commonly (and incorrectly) called : 'the failure of mass conservation', but, as 
stated before, m\ is more simply the probability for a vertex to belong to some finite cluster. 

On the other hand, the infinite set of Smoluchowski equations [Q] : 

~ir = 2 ^.i c * c i — K s ,jC s Cj , (3) 

i+j=s j 

are the coupled non-linear differential equations in the variables c s , i.e., in the concentrations 
of clusters of mass s. The time t includes both diffusion and reaction times, and these 
equations suppose irreversibility of aggregation, i.e., the cluster fragmentation is not allowed. 
The coefficients Kij represent the probability of aggregation between two clusters of mass 
i and j per unit of time. Some of them are explicitely known for different experimental 
conditions f|. But all such known aggregation kernels have the remarkable homogeneity 
feature : 

K = a x K 

for any positive a, with A called the homogeneity index. Maybe the simplest example of 
the homogeneous kernel is Kij = (ij)^. It has been shown theoretically that if fi is larger 
than 1/2, then there exists a finite time, say t cr , for which m' x becomes smaller than 1 for 
t > t cr ||. This can be interpreted as the appearance of an infinite cluster at the finite 
time t cr , and it is tempting to put in parallel the occupation probability p in the percolation 
model and the time t in the Smoluchowski approach. But we will see later that even if 
this parallel seems reasonable (p being the advancement of the aggregation process), some 
physical quantities behave quite differently. Note at last that in eq. @, as written above, 
the sum over j does not include the gel j = oo if any (since Cqo = l/oo = ), so the reaction 
between sol and gel is not taken into account. In principle, an additional term realizing the 
sol-gel aggregation should be added J7| on the right-hand member of eq. (|3|). 

II. REVERSIBLE AND IRREVERSIBLE SOL-GEL TRANSITIONS 

To see differences between both types of models, we shall treat two particularly simple 
cases : the bond-percolation on the Bethe-lattice, and the aggregation kernel K iy j = ij in 
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the Smoluchowski approach. They are chosen to be as close as possible in the sense that 
clusters generated by the percolation model are branched structures (without any loop) in 
an infinite-dimensional space. So, all their constituents are at the surface and reactivity 
must then be proportional to their mass |7|]. If the diffusion of clusters is negligible, e.g., 
when all clusters diffuse with the same velocity, or for highly concentrated systems, the 
corresponding reactivity kernels between two such clusters of mass i and j should be 
proportional to ij. 

The bond-percolation on the Bethe-lattice with coordination number z, has been solved 
by Fisher and Essam Here, the main result we are interested in, is the average concen- 
tration : 

c = z ^ Z ~ v s ~ l (I - v) {z ~ 2)s+z (A) 

° s Z ((z-2)s + 2)\s\ P 1 P) ' [) 

and the first normalized moment : 

with p* being the smallest solution of equation : 

p*(i- P *y- 2 = P (i- P y- 2 . (6) 

Let us define p cr = l/(z — 1). For p < p cr , the only solution of the above equation is : 
p* = p, but when p is larger than p cr , then there is a smaller non-trivial solution which 
behaves as p cr — \p — p cr \ near p cr . Above this threshold, the moment vn! x is smaller than 1 
and behaves approximately as 1 — 2(p — p cr ) /(l —p cr )- The marginal case z = 2 corresponds 
to the linear-chain case. 

Coming back to the concentrations, we can see that for large values of the size s, the 
following Stirling approximation holds : 

c s ~ s~ 5//2 exp(— as) , (7) 

with a given by : 



a = In ^(-^r 2 . (8) 




For this model, a power-law behaviour of the concentrations is seen at the threshold p cr , 
where more precisely c s ~ s~ T with r = 5/2. Outside of this threshold, an exponential 
cut-off is always present M . This sort of critical behaviour at equilibrium is analogous 
to the thermal critical phenomena, and in particular, there exist two independent critical 
exponents, for example r and a which is the exponent of the mean cluster-size divergence 
(here r = 5/2 and a — 1), to describe completely the critical features. 

The case of the Smoluchowski equations is quite different. Putting = ij in these 
equations, Leyvraz and Tschudi showed that there exists a critical value of the time, say 
t cr (here : t cr = 1), such that the solution is : 
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C s = 
C s = 

for size-distribution and : 



g s-2 

— —t 3 ' 1 exp(— st) for t <t, 
si 

— — exp(— s)/t for t > t cr 

m\ = 1 for t < t cr 

m'i = 1/t for t > t cr , 



(9) 



(10) 



for the first normalized moment. 

As explained in the previous section, the behaviour of w! x characterizes the sol-gel tran- 
sition, but some other features are interesting to compare to the percolation model. Firstly, 
one can see that the power-law behaviour is present for t > t cr and not only at the threshold, 
since for large s, we have : 

c s ~ s -5 / 2 exp(— as) ior t < t cr 



c s ~ s~ 5/2 /t for t > t cr , (11) 

with a = t — t cr — ln(t/t cr ). The whole distribution for the finite-size clusters evolves similarly 
and the appearance of a power-law behaviour is not a sign of the transition but rather a 
characteristics of the gelation phase. Secondly, it has been proved that for more general 
homogeneous kernels : Kij = (ij)^ , there exists a relation between the exponent r of the 
power-law behaviour and the exponent a of the divergence of the mean size, more precisely 
: t = a + 2. Here for = 1 we have r = 5/2 and a — 1/2. So, just one exponent is needed 
to describe the complete critical behaviour. In this sense, the reversible and irreversible 
sol-gel transitions, though close, are not equivalent. 



III. THE ORIGIN OF FLUCTUATIONS IN REVERSIBLE AND IRREVERSIBLE 

SOL-GEL MODELS 



As noticed by Einstein |Ipj , fluctuations of a macroscopic variable M at the equilibrium 
resembles a Brownian motion in the space of this variable and, hence, one expects the 
fluctuations to verify : 

<(M-<M>) 2 > ^ 

— i ~ Constant . (12) 

< M > K J 

'Constant' means here : independent of the mass of the system. This approach should be 
true for any short-ranged correlations, i.e., far from any critical behaviour. When close to 
the critical point, fluctuations are correlated throughout the whole system and Ornstein - 
Zernike argument yields [11] : 

< (M— < M >) 2 > ~ e »( d +n~ 2 ) , (13) 
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where e is the distance to the critical point, d is the dimensionality of the system, and v , 
7] are the two critical exponents related respectively to the divergence of the correlation 
length and to the divergence of the correlation function. In the same way, the average value 
of M behaves as < M >~ e 13 , so that, without making any additional assumption about 
relations between critical exponents, one obtains : 

<(M-<M>) 2 > ^ 

— i ——r 1 ~ Constant , (14) 

<M> S 

if and only if 5 = v(d + i] — 2)/ (3 . 5 is here a free parameter and should not be confused 



with yet another critical exponent. 'Constant' in ([LJ) means : independent of e for an 
infinite system. On the other hand, but by trivial finite-size scaling, this must also be a 
constant independent of the mass N of the system at the critical point. Moreover, if the 
standard relations between critical exponents hold, we obtain : 5 = 2 for any dimension 
d (for the Landau- Ginzburg theory we have to replace d by the upper critical dimension 
d c = 4 , above which the mean- field is valid). So the fluctuations of an order parameter 
in the thermodynamic systems are expected to behave differently at the critical point and 
outside of it. 

The case of irreversible aggregation models is quite different. Fluctuations in off- 
equilibrium physical processes are hard to analyze because they develop dynamically and, 
moreover, they often keep memory of a history of the process. In this sense, large fluctu- 
ations in such processes cannot be an unambiguous signature of the critical behaviour. In 
theoretical studies, it is often assumed that the fluctuations are irrelevant for the correct 
description of the mean properties of the system (e.g. in the Smoluchowski approach). This 
point will be revisited below. 

We will focus on the cluster-mass distribution n s , the cluster-multiplicity distribution 
P(Mq), and the distribution of the gel- mass for finite systems. The cluster- multiplicity M' Q 
is just the total number of clusters minus 1, i.e., the largest cluster is omitted. We wish 
to answer the following questions : (i) what is the importance of the power-law behaviour 
of the mass-distribution in detecting the critical behaviour? and, (ii) what are the scaling 
properties of multiplicity distribution and gel-mass distribution near the sol-gel transition 
for reversible and irreversible models ? The importance of the pertinent quantity M' Q in 
this context is that it is directly observable in many experimental situations where the 
system is not too large and the cluster masses are not directly accessible like , e.g., in the 
hadronization process in strong interaction physics or in the process of formation of the large 
scale structures in the Universe. Moreover, informations about the normalized moments m' 
or their fluctuations can also be obtained for large systems, such as in the polymerization, 
the colloid aggregation, or the aerosol coalsecence. 



IV. SCALING OF THE MULTIPLICITY-DISTRIBUTIONS AT THE 
REVERSIBLE AND IRREVERSIBLE SOL-GEL TRANSITION 

The multiplicity distribution is intensely studied in the strong interaction physics where 
simple behaviour of much of the data on hadron-multiplicity distribution seems to point to 
some universality independent of the particular dynamical process. Some time ago, Koba, 
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Nielsen and Olesen |L2| suggested an asymptotic scaling of this multiplicity probability 



distribution in strong interaction physics : 

< M> > P{M>) = *(*) , z = ^'m^* » ( 15 ) 

where the asymptotic behaviour is defined as < M' >— » oo, M' — * oo for a fixed 
M' / < Mq > ratio. < Mq > is the multiplicity averaged over an ensemble of independent 
events. The KNO scaling means that data for different energies (hence differing < Mq >) 
should fall on the same curve when < Mq > P(M' Q ) is plotted against the scaled variable 
Mq/ < Mq >. Extending this assumption, we suppose the more general scaling form : 

< M> > s P{M>) = *{z,) , z s = > > ( 16 ) 

with 5 a real parameter and $ a positive function. This form will be called the ^-scaling. 
KNO-case corresponds to 5 = 1. The normalization of the probability distribution P(Mq) 
and definition of the average value of Mq, provides the two constraints : 



/oo 
&(u)du = 1 
-<M£>i-« 

/oo 
u$(u)du = , (17) 
-<Af£>i-« 

which imply : S < 1 since $ is positive. 



As shown by Botet et al [T3[], the multiplicity distribution for the 3d-bond percolation 
model on the cubic lattice at the infinite-network percolation threshold exhibits the 5 - 
scaling with 5=1/2. Even though the system experiences the second-order critical phe- 
nomenon, fluctuations of the multiplicity-distribution remain small and the KNO-scaling 
does not hold. Of course, m' is not in this case an order parameter since r > 2 even though 
there is some irregularity in its behavior passing the threshold. This non-analyticity can be 
illustrated by the exact result for bond-percolation on the Bethe lattice. In this mean-field 
case, the normalized O tft -moment is : 

with : e = p — p cr , and e < 1. It is easy to see that there is a jump of the first p-derivative 
of m' : —z/2 for p — > p~ r , and (4 — 3z)/2 for p — > p+. The proper order parameter for this 
model is the normalized mass of the gel-phase, i.e., the mass of the largest cluster divided by 
the total mass of the system S max /N . Different probability distributions P(S max /N) for 
different system sizes N can be all compressed into a unique characteristic function (see 
Fig. 1) : 

< S max /N > P(S max /N) = n f S ^- < S ™ > ) , (19) 



which is an analogue of the the KNO-scaling function flU) ) of the multiplicity probability 



distributions. This new result is important since it seems to be a characteristic critical 
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behaviour of the order parameter. Now we will discuss what happens for a dynamical 
transition. 

We have simulated Smoluchowski approach by a standard Monte-Carlo binary aggrega- 
tion |14] . At each step of the event- cascade, a couple of cluster, with masses i and j, is 



chosen randomly with probability proportional to the time is then increased by the 

inverse of this probability, and the couple of clusters is replaced by a unique cluster of mass 
i + j. The largest cluster is always taken into account in this scenario. Let us discuss here 
the (fi = l)-case (K^j = ij). With the proper normalization, the critical gelation time is 
t cr = 1/N , where N is the total mass of the system. The power- law size distribution is 
indeed recovered numerically with the right exponent r = 5/2 (see Fig. 2). In contrast 
to the percolation model, the multiplicity distribution shows (5 = 0.2)-scaling (see Fig. 3) 
corresponding to very small correlated fluctuations. The scaling function, which in this case 
is asymmetric and sharp, is also quite different from the one for the percolation |13 . 



From the point of view of the criticality, the moment M[ should be a better candidate. 
The M{-distribution in the z$ - variable with 5 = 0.67 for masses iV = 1024, N = 4096 
and N = 16384 is shown in Fig. 4. This non-trivial value of 5 is here a signature of the 
transition : it disappears for t > t cr (see Fig. 5). But this is not yet analogous to the 
KNO-scaling because M[ is not exactly the order parameter. The true quantity to use is the 
reduced average mass of the largest cluster : S max /N = (N — M[)/N — 1 — m [ for which 
(5 = l)-scaling holds (see Fig. 6), as suggested by the phase-transition arguments, and in 
accordance also with the equilibrium percolation model. 

The results we have obtained here for the reversible or irreversible sol-gel transition 
in thermodynamical or dynamical systems are expected to satisfy the following general 
conjecture : 

• The occurence of (5 = 1)- scaling in the probability distribution P(M) of a cer- 
tain macroscopic quantity M is the sign of a critical behavior, with M as the order 
parameter; 

• The occurence of scaling with 1/2 < S < 1 in P(M) is the sign of a critical behavior 
in the system with M related to but not identical to the order parameter; 

• The occurence of 5 = 1/2 - scaling in the distribution P(M) is the sign that the 
variable M is not singular. 



These conjectures has also recently been checked for kinetic fragmentation systems ||15|| , and 
confort them. 



V. CONCLUSIONS 

We have studied here the scaling features of certain macroscopic quantities which are rel- 
evant for the sol-gel transition, such as the two first moments of the cluster-mass-distribution 
and the mass of the gel, i.e., the order parameter. We have employed two different mod- 
els : the percolation model which exhibits the equilibrium phase transition and the off- 
equilibrium Smoluchowski theory which shows a dynamical phase transition. The phase 
transition threshold in both cases is associated with the very particular scaling (|19"D of the 
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order parameter fluctuations. For Smoluchowski's equations, this scaling is seen uniquely at 
the critical gelation time t cr and disappears both for t < t cr and t > t cr . It is interesting 
to notice that for t > t cr , the cluster-mass-distribution is a power-law but , nevertheless, 
the order parameter fluctuations exhibit the small amplitude limit of fluctuations, i.e., the 
Gaussian scaling. 

For quantities other than the order parameter, similar scaling laws (see eq. (16)) are 
present but now with the different value of the parameter 5 of the scaling law. We have 
shown that the scaling : < M > s P(M) = $((M- < M >)/ < M > 5 ) for M,< M >^ 
oo, holds for any macroscopic quantity and the value of 5 is : 1/2 for the non-critical 
quantities, 0<5<1 (5^1/2 and 1) for critical quantities which related but which are 
non identical to the order parameter, and 1 for the order parameter, respectively. These 
new scaling laws could be used in the phenomenological applications either to find the phase 
transition threshold or to investigate whether the observed quantity is : regular, critical 
or whether it corresponds to the order parameter in the studied process. For example, in 
high-energy lepton-lepton collisions, the scaling behaviour with 5=1 has been reported 
in the multiplicity distributions of the produced particles for various collision energies [16| 
. This bahaviour has been interpreted as a signature of the asyptotic limit (the Feynman 
scaling) in the S'-matrix bahaviour, in spite of the fact that the Feynaman scaling has been 
experimentally disproved by the continuing rise of particle density in the central region. The 
analysis presented in this work offers an alternative and much more plausible explanation 
for these findings in terms of the off-equilibrium phase transition. It is indeed challenging 
to see the a priori unexpected relation between such different physical phenomena as the 
sol-gel transtition phenomenon and the particle production in the ultrarelativistic collisions 
of leptons. It would also be extremely interesting to apply the mathematical tools proposed 
in this work in the more typical aggregation phenomena as found in .... 
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Figure captions 



Fig. 1 

The gel-mass distribution for the 3d-bond percolation model on cubic lattices of different 
sizes : N = 8 3 (crosses), N = 10 3 (asterisks), N = 12 3 (triangles), and for the fixed bond 
probability p cr = 0.2488. Each point is the average over 200000 events. 

Fig. 2 

Double-logarithmic plot of the average mass-distribution of the Monte-Carlo Smoluchowski 
simulations with aggregation kernel = ij, for the two different system masses : 1024 
(diamonds) and N = 4096 (filled circles) , at the critical gelation-time (t cr = 1/iV). 

Fig. 3 

The multiplicity distribution for the Monte-Carlo Smoluchowski simulations with aggrega- 
tion kernel K it j = ij, for the two different system masses : 1024 (diamonds) and N = 4096 
(filled circles), in the ^-variable ( 5 = 0.2) , at the critical gelation-time. Each point is the 
average over 10 5 independent events. 

Fig. 4 

The M{-distribution for the Monte-Carlo Smoluchowski simulations with aggregation kernel 
— U> f° r the two different system masses : 1024 (diamonds) and N = 4096 (filled 
circles), in the zs- variable ( 8 = 0.67) , at the critical gelation-time. Each point is average 
over 10 5 independent events. 

Fig. 5 

The same as Fig. 4, but in the (5 = l/2)-variable and after the critical gelation-time 
(t cr = 2 /N). 

Fig. 6 

The ^m^-distribution for the Monte-Carlo Smoluchowski simulations with aggregation ker- 
nel K it j = ij, for the two different system masses : 1024 (diamonds) and N = 4096 (filled 
circles), in the KNO- variable at the critical gelation-time. Each point is the average over 
250000 independent events. 
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